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ABSTRACT 

We present a Bayesian technique based on a maximum entropy method to reconstruct the 
dark energy equation of state w{z) in a non-parametric way. This MaxEnt technique allows 
to incorporate relevant prior information while adjusting the degree of smoothing of the re- 
construction in response to the structure present in the data. 

After demonstrating the method on synthetic data, we apply it to current cosmological 
data, separately analysing type la supernovae measurement from the HST/GOODS program 
and the first year Supernovae Legacy Survey (SNLS), complemented by cosmic microwave 
background and baryonic acoustic oscillations data. We find that the SNLS data are compati- 
ble with w{z) = — 1 at all redshifts ^ z < 1100, with errorbars of order 20% for the most 
constraining choice of priors. The HST/GOODS data exhibit a slight (about Icr significance) 
preference for w > — 1 at z ~ 0.5 and a drift towards w > — 1 at larger redshifts, which 
however is not robust with respect to changes in our prior specifications. We employ both a 
constant equation of state prior model and a slowly varying w{z) and find that our conclusions 
are only mildly dependent on this choice at high redshifts. 

Our method highlights the danger of employing parametric fits for the unknown equation 
of state, that can potentially miss or underestimate real structure in the data. 
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1 INTRODUCTION 

' With the confirmation of the accelerated ex p ansion of the universe 
■ ("Perlmutter et al.' '1999"; lange et al.' '200 1'; 'Hoekstra et al.' '2002"; 
,Riess et al. 2004; Cole et al. 2005; Astier et al. 2006; Spergel et al. 

• [20061 ; iRiess et alJI2006h comes the inference that the cosmic dy- 
, namics are today dominated by a component that competes against 

• gravi tational collapse of matter and thus must hav e negative pres- 
. sure jprieman et ^ 1 19951 ; [Peebles & Ratral l2003h : this has been 

dubbed dark energy. All observations are presently compatible with 
dark energy being in the form of Einstein's cosmological con- 
stant A, a new form of matter-energy with equation of state (EOS) 
w — p/p = —1. However, it has been shown that an equation 
of state which changes with redshift, w{z), can mimic a cosmo- 
logical constant and fit the curre nt data if the parameterization of 
w is assumed to be a constant jLindejHooi ; Isimpson & Bridlel 
l2006l) . At the same time, an explanation based on the cosmologi- 
cal constant still suffers from the so-called "coincidence" and "fine 
tuning" problems, and it remains unclear whether selection effects 
of the kind embodied by anthropic arguments can offer a solution 
dStarkman & Trottall2006h . 

Determining whether dark energy is constant in time or has 
dynamical properties is one of the most pressing outstanding ques- 
tions in cosmology, as witnessed by the multiplication of ob- 
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servational proposals trying to elucidate the question (see e.g. 
iTrotta & Bowe 3i l2006h for an overview). We are thus led to ques- 
tion the simple solution of a constant w{z), especially in view of 
the fact that several recent works have highlighted the dangers of 
fitting current dat a by assuming a specific paramete rization for the 
dark energy EOS jBassett et aU2004lLindel2004) . 

The purpose of this paper is to explore the use of Bayesian 
statistical techniques based on a maximum entropy method to in- 
vestigate the time dependence of dark energy by imposing min- 
imal assumptions on the functional form of wl^z). Whenever ex- 
ternal (prior) information is used, its impact is clearly expressed 
by our formalism, making the reconstruction totally transparent. 
To do so we require information about the expansion of the uni- 
verse. The quality and quantity of observational data of cosmo- 
logical r elevance is rapidly incre asing: type la supernova (SNIa) 
(see e.g. lRiess et al.1 j2004r) ; lAstier _et al. (2006)) can be calibrated 
to serve as standard candles ("see lNugent et alj ((2006) for a recent 
proposal of using type II-P supernovae instead). Supernovae type la 
observations can thus measure the luminosity distance as a function 
of redshift, Dl{z), 

where the present-day Hubble constant is Ho = lOO/i km/s/Mpc, 
c is the speed of light in km/s and the redshift dependent Hubble 
function H{z) can be expressed in terms of the present-day matter- 
energy content of the Universe as (here and in the following we 
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assume a spatially flat Universe) 



+ Q.DE exp 



3/ ^-±^dz 
1 + z 



(2) 



Here, flm is the the density of matter in units of the critical density 
today and Qde is the (present-day) dark energy density in units 
of the critical density. The dark energy density time evolution is 
determined by its equation of state (EOS) 



Pde(z) = Pde{0) exp 



1 + w(z] 



dz 



(3) 



Observations of the luminosity-distance find a powerful geomet- 
ric complement in the use of "standard rulers" such as the position 
of the acoustic pe aks in the cosmic rn i crowave background (CMB) 
power spectrum ( iBennett et ai]|2003l : ISpergel et alJIZOOd) and the 
(transversal) baryon acoustic oscillation (BAO) si gnature recently 
measured in the galaxy matter power spectrum teisenstein et alj 
(2005); Cole et al. (2005)). Such data can be used to constrain the 
angular diameter distance Da (z) 



Da{z) = {1 + zYDl{z). 



(4) 



Apart from its geometrical impact on the angular diameter and 
luminosity distance relations, the properties of dark energy also in- 
fluence the growth of structures and can therefore be constrained 
through we a k lens ing (see e.g. Hu (2002); Hoekstra et al. (2005); 
Ijarvis et al ] ( l2006h ) and cluster counts data. In order to constrain 
w{z) from measurements of either Dl{z) or Da{z), we need to 
perform two derivatives, as it is evident from Q and ([5}. This is 
problematic if we consider the increase in the noise that accom- 
panies each derivative. In addition to information loss through this 
indirect determination of w(z), the current data tends to be sparse 
with a large sample variance. 

For these reasons, it seems timely and relevant to shift atten- 
tion to establishing more powerful statistical methods to extract in 
the most efficient and faithful way the information on w{z) con- 
tained in present and upcoming large data sets. The development 
of a technique that can cope with the patchy distribution in red- 
shift space while making minimal assumptions on the time prop- 
erties of dark energy is the logical next step towards improving 
our understanding of dark energy. In this paper we apply a modi- 
fied Maximum Entropy (MaxEnt) technique that has been diversely 
used to successfully reconstruct images and spectra under un- 
favoura ble conditions (for applications to astrophysical problems, 
see e.g . lBridle et alJ h998h ; lMarshall et alJ j2002h ;lMaisinger et al. 
( |2004) ). With a firm basis in probability theory, the method can 
be tailored to the needs of dark energy reconstruction from present 
data. Our application of MaxEnt aims at reconstructing the dark en- 
ergy EOS while minimizing our assumptions regarding the form of 
w(z). 

This paper is organized as follows. We firstly outline the sta- 
tistical framework of our technique in section |2] We then proceed 
to test our reconstruction method on synthetic data in section[3]and 
then apply it to present cosmological data in section!?] We offer our 
conclusions in section |5] 



MAXIMUM ENTROPY RECONSTRUCTION 
TECHNIQUE 



2.1 Motivation 

When attempting to constrain the nature of dark energy, a pro- 
cedure common in the literature is to Taylor expand the quantity 
Pde{z) or w{z) around z = 0, then constraining the expansion 
coefficients through the data ('Efstat hioi]|l999l ; iHuterer & Turnej 
^001; Weller & Albrecht 2002). An example of such a parameter- 
ization that is commonly emplo yed is w(z) = wq + (1 ~ a)wi 
where a 



(1+^) 



tChevallier & PolarskillioOll ; lLindej|2003l) . Al- 
ternatively, one might prefer to parameterize the time-dependence 
of the EOS using some smooth function (such as the ones used 
in lOick et al. ( 2006)), that is hoped will encapsulate the essen- 
tial features of the dynamics one wishes to constrain. Both pro- 
cedures are not free from the risk of giving misleading results, 
since they impose artificial assumptions on the form of the EOS, 
which often have no basis in any physical mechanism. Bass ett et al.l 
([2004) highlight the dangers of implementing such parameteriza- 
tions. Moreover, this will only be sensitive to departures from a 
constant density within a restricted set of models JPick et alJ200^ . 

iHuterer & StarkmanI f2003) introduced a principle compo- 
nent analysis (PC A) of the function w{z), with the aim of adopt- 
ing a parameteriz at ion appropriate to t he dat a sets used (see also 
iDick et al.1 l l2006h ; ISimpson & Bridld j2006h for related issues). 
These PCA modes are argued to form a natural basis in which to 
characterize dark energy evolution and by using only the first few 
well-determined eigenvectors in the reconstruction one tries to ex- 
clude noisy modes and thereby gain accuracy in the reconstruction. 
However, the method has the disadvantage of introducing an ill- 
controlled bias at high redshifts, i.e. the removal of strongly oscil- 
lating (and noisy) modes may mislead one to the conclusion that 
the EOS reverts to the fiducial model at large redshift with artifi- 
cially small error bars. While we recognize the merits of the PCA 
method, we wish to improve on it in this last respect by making 
fully explicit the assumptions that will control the behaviour of the 
reconstructed w{z) at large redshifts. 

The MaxEnt technique we employ has parallels with the well 
known maximum likelihood (ML) approach, but introduces new 
features ensuring that in the case where insufficient information 
is available the most likely distribution is the most uniform, i.e. 
the one with maximum entropy (or minimum information content. 
For an overview of t he connection b etween entr opy and informa- 
tion content, see e.g. iTrottal ( |2005|) ; iKunz et alj fcood) ). In other 
words, where ML merely maximizes the likelihood, often unneces- 
sarily overfitting the noise in the data, MaxEnt seeks the optimum 
trade-off between a smooth, maximally entropic distribution and 
the rough distribution mapped out by the data. The most impor- 
tant characteristic is that the MaxEnt method is auto-regulating, 
i.e. the amount of smoothness (or raggedness) in the reconstruction 
is consistently determined through the data themselves (see section 
l2.2.4l below). In our Bayesian perspective, extra information com- 
ing from prior beliefs or theoretical prejudice can be naturally in- 
corporated in the reconstruction via the MaxEnt prior. As we show 
below, this gives MaxEnt the power to cope with situations where 
the dimensionality of the parameter space potentially exceeds the 
number of data points, a difficult reconstruction problem that is ill- 
defined under ML techniques. This feature clearly makes MaxEnt 
highly applicable to the case of dark energy reconstruction. 
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2.2 The MaxEnt formalism 

The task at hand is to determine the EOS of dark energy from sparse 
data on Dl, Da and, to a limited extent, on H{z) (as encapsulated 
by today's detection of BAO). We consider a piece-wise constant 
w{z) in A*' bins out to a maximum redshift Zmax- Let Wj be the value 
of the EOS in the 7-th bin, 1 ^ j ^ A'^. In analogy with the treat- 
ment given in lSkillind ( 1 19891) for the case of image reconstruction, 
we gather all the EOS bins values in an "image vector" w. We seek 
to determine the posterior distribution of w given the observed data 
D, Pr{w I D). This is obtained through Bayes' theorem as 

Pr(D I w)Pr(w) 



Pr{w I D) 



(5) 



Pr-(D) 

The quantity Pr-(w) is the prior probability representing all the in- 
formation about the distribution w before the data D has been col- 
lected; Pr(D I w) is the likelihood and describes the underlying 
statistical process and Pr(T)) is the model likelihood ("evidence"), 
which is relevant for model selection questions but that is unim- 
portant in this case. We shall therefore neglect this proportionality 
constant in the following. 

2.2.1 The MaxEnt prior 

The principle of MaxEnt is employed to determine a prior Pr (w) 
that encapsulates all the external information abo ut w{z) we wis h 
to specify in the absence of the data. Following ISkillind dH), 
we adopt the principle that the least biased model that encodes any 
given prior information is the one which maximizes the entropy 
of the distribution while remaining consistent with the informa- 
tion. This prior is appealing for its characteristic of maximising 
the uncertainty (entropy) of the distribution th us making minimal 
assumptions. The MaxEnt prior i lSkillindll989h takes the form 

exp (aS'(w, m)) 



Pr (w I a, m) — 



Zs 



(6) 



where S (w, m) is the entropy of w relative to the model m and a 
is a regularizing constant. The model m defines the image vector 
to which w reverts in the absence of any data, and as such it defines 
a measure on the DE parameter space. The information entropy is 
analogous to the thermodynamic entropy in statistical mechanics, 
which is given by the logarithm of the number of states by which 
one can arrive at a given macroscopic constraint. In the same way, 
the information entropy can be described as the logarithm of the 
numbers of ways in which one can arrive at a pa rticular w in a Pois - 
son process when starting with the model m dSilver et alj[l990h . 
The entropy S for an A'^-dimensional discrete parameter space is 
(ISkillinjl989[) 



S(w,m)=^ 



Wj log 



(7) 



The log term is reminiscent of the Kullback-Leibler divergence be- 
tween w and m, encoding the amount of information present in w 
with respect to the model m. In our case, we do not apply the en- 
tropy to the values of w directly, but rather to the space of coeffi- 
cients of an expansion of w in a series of basis functions (that we 
choose to be top-hat functions in redshift space, see section [231 for 
details). We can think of the coefficients of the expansion as a series 
of weights that encode how much each basis function contributes to 
the total wi^z). We can then apply the MaxEnt prior on the space of 
these weights, by thinking of them as expressing the relative contri- 
bution of each basis function - in other words, in a phenomenolog- 
ical way we take the weights to represent relative probabilities for 



the presence of each basis function in the final w{z) . Below, we will 
use the notation w as a shortcut to indicate the vector of weights of 
the dark energy expansion. The same applies to the model m, that 
in the entropy term is represented by its expansion coefficients in 
the chosen basis functions. 

Evidently, S'(w) (for a fixed m) is a convex function which 
reaches a maximum for w = m with a value 5* = 0. Thus in the 
absence of any information from the data, the entropy term reverts 
to the model. The normalizing partition function for the entropy is 
given by 



/ 



Zs ~ I exp aS{w) det[g]^^^d'^w 



(8) 



The measure is defined as the invariant volume det g^^^ of the met- 
ric g defined on the space where gu = 1/wi and gij = for i ^ j 
(also known as the Fisher information matrix). By expanding to 
second order around the model w = m (at the maximum S = 0), 
we obtain the partition function in the Laplace approximation: 



Zs = ( 



iV/2 

2tv) 



(9) 



where A'^ is the number of parameters, in our case the number of 
expansion coefficients for 

2.2.2 The likelihood 

The likelihood is defined as of the probability of the data given the 
parameters: 



Pr(D I w) = 



Zl 



(10) 



and is the probability that the observed data D could have been 
generated from a given w. For data D subject to Gaussian noise 
the likelihood function is 

£(w) = i(D - fMf[C-']{T> - /(w)), (11) 

where C is the data covariance matrix and /(w) denotes the func- 
tional dependence of the observable on the DE parameters in our 
case, f = H ov f = Da (Da and Dl being simply related through 
the redshift, see llU). In the case of independent data points with 
uncorrelated noise, the covariance matrix is diagonal with the non- 
zero elements being the variances of each measurement, denoted 
hy af, i — 1, . . . , No. The normalizing partition function for C is 



Zl 



exp(£(w))d^°D 



(12) 



Using the identity for the normalized probability distribution we 
obtain: 



{2tt 



Vdet[C7-i] 



(13) 



2.2.3 The posterior probability 

From the likelihood in Eq. dlOb and the prior in {6}, we obtain from 
Bayes' theorem the posterior probability for the DE parameters w: 



Pr{w I D, a, m) oc exp(aS'(w) — £(w)). 



(14) 



Given that £(w) is quadratic in w and ^(w) is a convex function, 
the above is well-constrained with the peak of the posterior for w 
being determined by a competition between S and C, mediated by 
the value of a. We thus see that the MaxEnt prior will be useful in 
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the case where the parameter space dimensionaUty exceeds the size 
of the dataset in that the entropy is incorporated as a regularization 
to avoid over-fitting, while retaining maximum flexibihty in the un- 
derlying parameterization of w(z). S penalizes the excess "struc- 
ture" in the data, with the regularizing parameter a dictating the 
degree of this smoothing. The choice of a is thus very important: a 
small value of a produces little smoothing and the likelihood will 
dominate, resulting in a reconstructed distribution where the noise 
might be mistaken as real structure. Alternatively, too large a value 
for a leads to information loss, with the entropic prior overriding 
the information coming from the data. 



2.2.4 The regularization parameter a 

In order to select the correct value of a we add it to the hypothesis 
space as an additional parameter and let the data select the optimum 
value. Using once more Bayes' theorem we have 



Pr {a j D, m) oc Pr (a) Pr(D | a, m), 



(15) 



where Pr{a) is the prior on a. The joint posterior probability is 
then (up to irrelevant constants) 

Pr(w,a,D|m) oc Pr{a)Pr{w \ a)Pr(D | w) 

Zs(a)ZL 

We adopt a Jeffreys' prior on a, which is flat in 7 = log a, re- 
flecting ignorance of the scale of the variable, within the range 
— 10 ^ 7 ^ 10. This corresponds to choosing Pr{a) oc In 
the final inference on w we marginalize over the nuisance param- 
eter a, even though the distribution of a is usually fairly strongly 
peaked and thus marginalization is almost equivalent to maximisa- 
tion (i.e., simply fixing a to the value of the peak of the posterior). 



2.2.5 Model specification 

The joint posterior in Eq. l ll6t is conditional on the specific choice 
of model m, to which the reconstructed w defaults in the absence 
of data. The entropic prior distribution is introduced to penalize the 
posterior for unwarranted complexity. Given that the model is the 
vector to which w should revert in the absence of data, it must rep- 
resent maximal smoothness. We need to establish what distribution 
m will encompass this in the context of the EOS of dark energy. 
In image reconstruction the default model is usually taken to be a 
flat surface equal to the mean of the data. When the data is then in- 
cluded via the likelihood, variation about this mean is introduced. 
In our case, this means choosing a contant model, m — const that is 
uniform in redshift space. There is however no obvious choice for 
the magnitude of this constant. There are various possible choice 
of m, reflecting different prior beliefs about the dark energy EOS. 
One can thus usefully think of m as encapsulating a fiducial, refer- 
ence model we want to test the observations against. One possibility 
is to set m = — 1 at all redshifts, thus representing a cosmological 
constant. This is recommended if we are testing for deviations from 
w{z) = —1: if significant deviations from the model are found in 
the reconstructed EOS, this is an indication that the data are in- 
formative enough to override the entropic pull towards the model, 
and thus that such deviations are likely to reflect real structure in 
the data. A more skeptical attitude towards dark energy might be 
encapsulated by choosing a constant model m = 0, which corre- 
spond to a pressureless, dust-like fluid. In this case, if the recon- 
structed 'w{z) assumes values below 0, this can be interpreted as a 



strong indication for the presence of a fluid with negative pressure, 
with data being strong enough to dominate the entropic tendency 
for a pure matter Universe. In principle, a theoretical prejudice in 
the form of a particular redshift-dependence of m could also be 
implemented easily in the same fashion. 

Finally, one can also employ Bayes theorem to take m into 
the joint posterior, by writing 



Pr{w, D, a, m) oc Pr{m)Pr{w, D, a|m) 



(17) 



and marginalizing over m in the left-hand-side, after specifying 
a prior over the model space, Pr(m). In this work we take all the 
model vectors to be constant over the whole redshift range, thus the 
specification of the model amounts to the choice of the value of the 
constant. We restrict our considerations to the range — 1 ^ m ^ 0, 
and when performing a marginalization over the model we will take 
a flat prior in this range, i.e. Pr(m) = const. 



2.3 Dark energy parameterization and reconstruction 

As motivated above, we decompose 'w{z) into a weighted sum of 
orthogonal functions in redshift space, with the parameters being 
given by the weights encoding the amount that each function con- 
tributes to the overall 'w{z). There are of course several different 
meaningful expansion functions such (for example, principal com- 
ponents, Chebychev functions, etc) but we make use of the sim- 
plest option, decomposing w{z)is into a series of A*' step-functions 



-2 + ^Ci$,(2 



(18) 



where $i{z) = 1 for Zi — Az/2 < z < Zi + Az/2 and $i{z) = 
everywhere else. Since the least stringent limits we will impose 
on w{z) are — 2 ^ to^ ^ 0, the above ensures that the expan- 
sion coefficient d (i — 1, . . . , N) are positive numbers, a neces- 
sary requirement for our entropic prior. The parameter space w is 
thus constructed from the co-efficients d themselves, which are 
allowed to vary within the range ^ Ci ^ 2. An advantage of 
this piece-wise parameterization of is that it will be possible 
to capture a sharp change in the EOS, provided the binning is suffi- 
cient. In order to capture different features of the time evolution of 
the EOS, other expansion functions may be more appropriate. We 
experimented with Chebychev functions and found that their oscil- 
latory behaviour was not helpful in reconstructing sharp changes in 
the EOS. Such smooth functions might be more useful if one wants 
to test quintessence models exhibiting a gentle evolution of w(z). 

In the bulk of recent analyses the limit — 1 w ^ is 
imposed; the lower limit stemming from the null energ y condi- 
tion which must be satisfied for dark energy to be stable jAlcanij 
2004). Although models of dark energy that allow w < —1 vi- 
olate the weak energy condition in the context of general relativ- 
ity, th ese "phantom" c omponents have been studied by many au- 
thors ( ICaldwellll2003) . There have been claims that such phan- 
tom behaviour is unstable when regarded as a quantum field the- 
ory (Carroll et al. 2003). From a phenomenological perspective, it 
makes sense to both restrict the range of our reconstruction to lie 
within — 1 ^ w ^ 0, but also to extend the parameter space to 
values below w = — 1 to check the stability of the reconstruction. 
We will thus presents results also for the case where the equation 
of state can attain values as low as w = — 2. 

In this work we assume flat spatial sections, and thus f2m ~ 
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[t] 



Parameter 



Prior 





0.0. . 


1.0 


Top-hat 


O,^ 







Flatness imposed 


Ci 


0. . 


2 


Top-hat 


7 = log(o) 


-10. 


. 10 


Top-hat 


m 


-1 . 


.0 


Top-hat 



Table 1. Priors on cosmological and nuisance parameters used in the anal- 
ysis. We employ a Jeffreys' prior on a, i.e. we take the prior to be flat in 
7 = log a to reflect ignorance on the scale of the regularizing parameter a. 



1 — Qde- The parameters included in the hypothesis space are sum- 
marized in Table [T] These are Ude, the Hubble parameter today, 
Ho in km/s/Mpc, and the coefficients of the DE decomposition, d, 
as described by Eq. ilSl . This generic characterization requires the 
number of expansion functions (which can be effectively character- 
ized as top-hat bins) to be sufficiently large for this to be a suitable 
description of As described above, we also include the nui- 

sance parameters 7 — log a and the value of the model EOS, m, 
whenever this is marginalized over. 

Assuming uncorrelated Gaussian noise, the log-likelihood of 
a point in parameter space is given by 



2iog/:(w)^^( ^'--^("'-) )', 



(19) 



where for each datum i at redshift Zi we have f{zi,w) = H{zi) 
for future radial baryonic oscillation measurements, f{zi,w) = 
Da (zi) for present and future transversal baryonic oscillation data 
and CMB data and f{zi,w) = DL{zi) for SNIa data. Further- 
more, is the measurement variance. The Hubble parameter as a 
function of redshift is obtained via Eq. |2}, where the energy den- 
sity is calculated using, for Za falling within the i-th bin; 

3(l+mi) 



pDEjZa) 

Pde{0) 



1 + Za 



n 



1 + Z^ ~ Az^/2^ 



J = l 



(20) 



1 + zj - Azi/2^ 

From the above, the angular diameter distance DA{za) can 
then be computed using equations illH . For piecewise constant 
w{z) we employ the trapezoid rule to approximate the integral, ob- 
taining 



DA{Za) = 



+ 



1 



1 + Za 
1 



Sz 

y 



1 



H{2 



(21) 



Sz 1 

a-l 

+ E- 

The binning of the integral, defined by Sz is determined based on 
a fixed level of fractional accuracy for the integration, that we set 
to lO"^'^ as determined by the extrapolation error estimate. Finally, 
the entropy of a vector w with respect to a model m is given by 



5(w) = Yl 



Ci - Mi ~ C\ log 



(22) 



where the Mi, i = 1, . . . , N are the coefficient of the model m in 
the expansion dlSt . 

To sample the posterior probability distribution efficiently we 
use a Monte Carlo Markov Chain (MCMC) which employs a 



Metro polis algorithm. Fo r more details about MCMC, see e.g. iNeall 
( I1992I) : [Lew is & Bridii j2002h . Since the MaxEnt method is de- 
signed to achieve the optimal reconstruction independently of the 
number of degrees of freedom in the parameterization of we 
expect that the number of basis functions A*' will not affect greatly 
the reconstruction, as long as A'^ is chosen large enough to cap- 
ture the possible structure in the data. In the following we choose 
TV = 10 but we have checked that the results do not vary much if 
one uses A'^ = 5 or = 20 instead. 

When using actual data, we divide the redshift range spanned 
by either the SNLS or the HST/GOODS supernova measurements 
into A'^ = 10 equally spaced bins, corresponding to the A' basis 
functions for wiz). We then extend the last bin to cover all of the 
redshift range to last scattering when computing the angular di- 
ameter distance to the CMB. In other words, we take w{z) to be 
constant (but not fixed to —1) between the redshift of the high- 
est supernova in the samples and 2 — 1089. This extr apolation is 
weake r than the 'strong' prior used in the analysis of iRiess et al.l 
( l2006h . which assumed that w = —1 at 2 > 1.8. 



3 DEMONSTRATION OF THE MAXENT METHOD 

We now proceed to test our MaxEnt reconstruction method with 
synthetic data. Our benchmark dataset consists of Nd = 10 mea- 
surements of H{z) and Nd ~ 10 of Da{z) (or equivalently, 
Dl(z)) distributed uniformly in the redshift range ^ 2 ^ 1. 
Although the actual measurements will in reality be less homo- 
geneous, this does not represent a problem for our reconstruction 
algorithm, as we show below. Existing measurements of Dl{z) 
out to 2 ~ 2 will be vastly improved when future surveys such as 
DES or LSST will be able to observe thousands of SNIa per year 
( Abb ott et al . 2005; Tvson 2006) and space-based projects such as 
SNAP jAldering et alj2004) . ADEPT or DUNE will provide obser- 
vations at redshift beyond 2 ~ 0.8. Future spectrographic surveys 
such as the Wide-Field — Multi-Object Spectrograph (WFMOS) 
or HETDEX ought to deliver constraints on Da (2) of 1% a t 2 ~ 1 
and 15% at 2 ~ 3 (la ) and H{z ~ 1) to 1.2% iGlaze brook et al.1 
l2005l ; lKelz et al.l2006h . with better performance still to be expected 
when the Square Kilometer Array will come online jBlake et"ail 
[2004). 

We add Gaussian noise to our synthetic data as a fixed percent- 
age of the true value of the measured quantity. We use an optimistic 
noise level of 1%. Although it is beyond the scope of this paper to 
make quantitative predictions about the performance of future sur- 
veys in reconstructing the EOS, our benchmark data set roughly 
reflects the potentiality of future observations. We also test the per- 
formance of the method when the signal-to-noise level is degraded 
by a factor of 10, in order to check for bias in the reconstruction due 
to our entropic prior when the quality of the data is poor. In this case 
we use a noise level of 10% in the luminosity distance and Hubble 
function measurements. 

We show in Figure[T]the reconstructed EOS for our benchmark 
scenario with high quality observations (a = 1%, Nd ~ 20 ob- 
servations). In all three panels, we have marginalized over the prior 
model m. We notice that the reconstruction is satisfactory in all 
three cases. We have checked that the w extracted when marginal- 
izing over the prior model m has comparable accuracy to the case 
of a fixed model m. The bottom panel shows how the method deals 
with gaps in the redshift range of the observations: the smoothing 
effect of the entropic term enlarges the errors in the region where 
no data are available, while the reconstructed EOS tracks the true 
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Figure 1. Reconstructed EOS w (black errorbars showing la posterior con- 
straints) using our MaxEnt method for high-quality synthetic data {Nn = 
20, a = 1%, location shown by the blue diamonds) in the redshift range 
2 1. Top panel: the true EOS (blue, dotted line) is taken to be a 
step function. Middle panel: the true EOS shows a sharp peak at z ~ 0.4. 
Lower panel: the true EOS is slowly evolving with redshift, and the syn- 
thetic data are now clustered at low and high redshift. Despite the absence 
of data points at intermediate redshifts, the high-z reconstruction tracks the 
true function with reasonable accuracy, while the intermediate redshift er- 
rors increase con'espondingly. In all three cases, the value of the prior model 
m for 111(2) has been marginalized over and the MaxEnt reconstruction sat- 
isfactorily recovers the true EOS. We have plotted horizontal lines at ui = 
and = — 1 to guide the eye. 
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Figure 2. Reconstructed EOS w (eiTorbars showing Icr posterior con- 
straints) for noisy synthetic data {Njj = 20, a = 10%, location shown by 
the blue diamonds). Top panel: the reconstructed EOS for a model m = — 1 
(red errorbars) has collapsed towards the model due to the entropic prior for 
redshifts 2 ^ 0.3, while the reconstruction with m = (blue en'orbars) 
tracks better the true EOS (blue, dotted line). Bottom panel: after marginal- 
ization over the model value, the bias in the reconstruction has disappeared, 
but the en'orbars have become suitably larger. 



value at small and large redshifts, where the data are clustered. Here 
we have employed a 10-dimensional w, but we have checked that 
increasing the number of elements to A'^ = 20 does not lead to any 
significant change in the reconstruction. As expected, the error bars 
in the regions where 'w{z) is well constrained by observations are 
considerably smaller. 

We now turn to the case where the data are noisy, and hence 
we expect our entropic prior to play a more important role in the 
reconstruction. The top panel of Figure |2] shows the reconstructed 
EOS for a noisy data set of Nd ~ 20 measurements with noise 
(J — 10% and a slowly evolving true EOS, w{z) = 1 — ln(l + 2). 
In the top panel, we show the result when employing as prior model 
values the two liming cases of w with strong theoretical preju- 
dice: m = — 1 (cosmological constant, red errorbars) and m = 
(Einstein-de Sitter Universe, blue errorbars). In the high z bins the 
reconstruction becomes increasingly mismatched with the under- 
lying true EOS. Because the dependence of the data on the EOS 
requires integrating the EOS over redshift, any error in the recon- 
struction at low redshift is accumulated as z increases. As a result 
the entropy tends to dominate over the likelihood and the mean pa- 
rameter values collapses towards the model at higher redshift, espe- 
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cially for the case where m = — 1. Even though the reconstruction 
has appeared to degrade for the m = case, it is encouraging 
that the mean values of the parameters in the lower redshift bins 
{z < 0.3) are reasonably close to the true values. 

Evidently the choice of prior model m does have some bear- 
ing on the reconstructed value of the parameters at high redshift 
(given that the entropy dominates the posterior for poorly infor- 
mative data) and this must be kept in mind when interpreting the 
results when dealing with noisy data. This example highlights the 
problem of distinguishing a genuine affinity for a certain function 
that happens to closely resemble the model from a strong default 
towards the model on account of noisy observations, i.e. how will 
we interpret a result very close to w{z) = —1? In this case, there 
are two options: the observational evidence is noisy and uninforma- 
tive, leading to entropy domination, or the data is good and favours 
an actual value close to —1. In the latter case, attempting a recon- 
struction with m = will allow to test the strength of the data in 
pulling w{z) towards the cosmological constant value. An alterna- 
tive mean of recovering w{z) in a truly model-independent way is 
to include the elements of the model vector in the hypothesis space, 
as discussed in section [2.2.5l Marginalizing over m then amounts 
to selecting the optimum distribution at each sampled point. The 
model still represents the most entropic distribution given that it is 
uniform in z-space. The result of this procedure is shown in the 
bottom panel of Figure |2l where marginalization over the model 
has cured the skew in the reconstruction observed above, albeit at 
the price of delivering larger error bars. 



This becomes a standard ruler because the scale of these acous- 
tic os cillations is self-calibrated under standard recombination ( 
l2005h . It depends solely on the photon-baryon ratio and radiation- 
matter ratio at recombination which are determined with excel- 
lent precisio n in the CMB power sp ectrum from the CMB peak 
morphology jEisenstein & W hite 2004). The change in the appar- 
ent size of this scale from recombination to the present will de- 
pend on the expansion history of the universe through projection 
effects. These acoustic features appear as rings in angular and red- 
shift space ( Hu & Haiman 2003). The actual measurement from the 
SDSS LRG sample is of the dilation factor, defined as 



Dv{z) = 



Dl{z) 



H{z) 



1/3 



(24) 



where the comoving angular diameter distance Da is taken as the 
transverse dilation and the line of sight measurement of this scale 
Sf-\ is taken to be the radial dilation. The observed correlation 

H(z) 

scale constrains a function of the dilation factor, and a single data 
point is reported in lEisenstein et a Di l2005h : 



A = Dv-(0.35) 



H"2 



0.35c 



0.469 ±0.017. 



(25) 



(see also ICole et al.l bOOSh for a similar detection of the acoustic 
feature in the 2dF catalogue). The above assumes Vlbh^ = 0.027. 
The log-likelihood for the baryon acoustic oscillation data is given 
by 

{A ~ 0.469)2 



2 

Xbao 



0.0172 



(26) 



4 APPLICATION TO PRESENT-DAY DATA 

We now apply our reconstruction procedure to actual data, en- 
compassing cosmic microwave background observations, baryonic 
acoustic oscillations measurements, determination of the present 
value of the Hubble parameter and two different supernovae type 
la samples. 



4.1 Data sets 

4.1.1 Baryon Acoustic Oscillations 

The comoving sizes along the line of sight, r|| (in redshift space), 
and in the traverse directions, rx, of a feature sitting at a redshift 
z are related to the redshift range Az covered and the angle sub- 
tended AO, respectively, by 

cAz 

rii^-^, t^^{1 + z)Da{z)^Q. (23) 

If the absolute values of ri_ and r|| are known, they become stan- 
dard rulers giving us a handle on H{z) and Da {z) . If only the rela- 
tive sizes are known, then the standard rulers are expressed in units 
of Hq. If only the ratio ^n/rx is known, this becomes the Alcock- 
Paczynski test. The baryon acoustic oscillation phenomenon can be 
used as such a standard ruler 

After recombination, when the Universe becomes neutral and 
photons free-stream from the cosmic plasma, the driving force of 
the harmonic oscillation is removed and the sound speed of the 
now-neutral medium essentially falls off to zero, ending wave prop- 
agation. The spatial distribution of the baryons at this stage will 
then reflect the characteristic scale of the acoustic waves. Seeing 
as the perturbations in the baryon and CDM distributions seed the 
formation of large scale structure, we exp ect to see acoustic peak s 
in the late-time matter power spectrum l lEisenstein et alj l2005h . 



4.1.2 The Supernova (SNIa) data 

Type la supernovae are good candidates for standard candles and 
are useful in determining distances on extragalactic scales. Due to 
the complexity of the physics involved, the SNIa are not perfect 
standard can dles, having a disp ersion of 0.3 — 0.5 mag in their peak 
magnitudes ( Istraumaruil2006l) . However the peak brightnesses ap- 
pear to be tightly coiTelated to the time-scale of their brightening 
and fading and one can extract an empirical relation between ab- 
solute peak luminosity and the morphology of the light curves to 
constrain the absolute brightnesses, and thus obtain measurements 
of Dl{z). 

The first group of superno vae, termed the 'g old' set, is from 
the HST/GOODS programme l lRiess et alj|2004j) . complemented 
by the recently di scovered higher redshift supernovae, reported in 
iRiess et alj ( I2OO6I) , while the second sample is taken from the Su- 
perno vae Legacy Survey (SNLS) iAst ier et al. 2006). As discussed 
in e.g. lWang & Mukheried j2006l) , it appears that there are system- 
atic differences between these two data sets that arise from differ- 
ences in the data processing. It is therefore necessary to consider 
the two datasets separately, and compare the results as a consis- 
tency check. 
The 'gold' sample: 

The distance modulus /i is defined as the difference between the 
apparent magnitude m and the absolute magnitude M : 



fj, = m ~ M = 5 logj^ 



Dl{z) 
lOpc 



(27) 



Given that the absolute magnitude M is a unknown, we consider 
the following distance modulus iii{8M): 

fii{5M) ^ fif ~5M. (28) 

Here 5M is the difference between the mean of the true absolute 
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magnitudes and the estimated absolute magnitude, while the are 
the observed magnitudes after dust corrections and recalibration 
through the shape of the luminosity evolution function. The quan- 
tity 5M is the difference between the true mean absolute magni- 
tude and the estimated absolute magnitude of the supernovae and 
is marginalized over. 
The SNLS sample: 

The SNIa data from the Supernovae Legacy Survey (SNLS) are 
reduced in a different manner in that the light curves provide con- 
straints on various parameters which are then used to calculate the 
effective apparent m agnitude. For a description of the calculation 
of iii see lPick et alJ (.2006) . From the observed /li with variances 
(jjj i for each set of SNIa, we perform an analytical marginalization 
over the absolute magnitude M. Defining the quantities 



Mi 



1 



and 



N 



(Slogio DL(Zi) - 



(29) 



(30) 



(31) 



the M-independent log-likelihood for the SNIa is calculated as 

xliv = /i-^. (32) 
c 

4.1.3 CMB and HST data 

The WMAP 3-year measurement of the CMB shift parameter 
describing the location of peaks in the CMB power spectrum 
serves to constrain t he angular diameter distance to last scattering 
dSpergel et alj|200^ . This is independe nt of most assumptions of 
the form of dark energy, and is given bv dWang & Mukherieel2006h 



R = Q.U'^H, 



''f 

Jo 



dz 



1.70 ±0.03, 



(33) 



where zcmb is the redshift to last scattering, taken in our case to 
be 1089. 

We also include the constraint on the present value of the Hub- 
ble constant obtained by the Hubble Space Telescope Key Project 
iFreedman et al by using a Gaussian likelihood with mean 

and standard deviation given by Ho = 72 ± 8 km Mpc s~^. 



4.2 Results 

We plot in Figure [5] the results of our reconstruction from the 
SN type la data from the SNLS and from the HST/GOODS pro- 
grammes. In both cases we have added the CMB, HST and BAO 
measurements, and we have marginalized over the model, in order 
to be as conservative as possible. Furthermore, the maximum range 
for the EOS has been taken to be —2 ^ w ^ 0. We plot regions 
encompassing 68% of posterior probability for each w bin - notice 
that since these are marginalized values, their magnitude is inde- 
pendent of the correlations between reconstructed points (the issue 
of correlations is addressed in detail below, see Figure©. 

We do not find any significant deviation from a cosmological 
constant behaviour from the SNLS data (see top panel of Figure[3]i 
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Figure 3. Reconstructed EOS (marginalized errorbars encompassing 68% 
of posterior probability) from SNLS data (top panel) and HST/GOODS data 
(bottom panel), including CMB, HST and BAO measurements, as well (no- 
tice that the redshift range is different for the two panels). The prior model 
m is constant in redshift and has been marginalized and the assumed range 
of the EOS is — 2 ^ w ^ (both are conservative choices). The horizontal 
lines indicate the upper bounds of the allowed w range (solid) and position 
of lii = — 1 (dotted) in order to guide the eye. The SNLS data do not show 
significant deviations from w = —1 over the whole redshift range. The 
reconstruction from the HST/GOODS data is also consistent with a cos- 
mological constant, although it appears to slightly prefer a higher value at 
redshift z ~ 0.5. The best-fitting points are all very close tow = —1. 



at all redshifts. Our posterior constraints in this case are rather con- 
servative, as a consequence of the assumptions made in the recon- 
struction (i.e., large w range and marginalization over the model). 
When using the HST/GOODS data, the recovered w in the first bin 
agrees with that found for the SNLS data. Given that a number of 
the SNIa in this bin are common to both surveys, this provides a 
consistency check. The reconstructed EOS from the HST/GOODS 
data is however found to prefer slightly higher values in the third 
bin {z ~ 0.5), excluding the cosmological constant value to a 
little bit more than la significance. The significance of this rise 
however has to be assessed with care, especially if we recall that 
the constraining power o f the SNIa data degrades in this region 
jSimpson & Bridlell2006h . At redshifts above z ^ 0.5, the mean 
of the reconstruction settles around w ~ —0.7, although we notice 
that the best fitting points remain very close to «; = — 1 (red trian- 
gles in Figure[3}. The implied early-time behaviour of dark energy 
is con sistent with the result of «; = — 0.8l']' o found in lRiess et al.l 
( l2006i) for z > 1, using what they call their 'strong' prior. If instead 
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Figure 4. Correlation matrix for tlie reconstructed EOS parameters from 
HST/GOODS (upper left) and SNLS (lower right) data sets, including 
CMB, HST and BAO measurements. The strongest (anti)correlations are 
between values in bins 2 and 3 which roughly coincide with the redshift 
position of the BAO measurement. 



the EOS is assumed to be constant over the entire redshift space 
(i.e., if we reduce the number of w components to A'^ = 1), then we 
obtain from the HST/GOODS d ata w = -0 .89 ± 0.07, in agree- 
ment with usual results (see e.g. iRiess et al ] ( 120061 ). This clearly 
demonstrates the danger of assuming 'w{z) to be time-independent, 
as one would miss in this way possible features in the data. 

The use of the entropic prior introduces c orrelat ions among 
the reconstructed points (see lHuterer & Cooravl dlOOSh for a tech- 
nique to extract uncorrelated band power estimates of the EOS). 
The correlation coefficients from the posterior distribution over the 
w parameters are shown for both data sets in Figure |4] We notice 
that correlations are in general relatively mild, flattening around the 
level of ~ 20% for correlations with bins at larger redshifts, where 
the entropic prior becomes more important. The strongest correla- 
tions (at the level of ~ 50%) are observed among parameters in the 
2nd and 3rd redshift bins, where the BAO measurement strongly 
constrains the EOS and due to the fact that the observables are in- 
tegrated over redshift, we expect a negative correlation among the 
well-constrained value at the position of the BAO measurement 
and the w values at lower redshift. 

We now investigate the case where we impose that w ^ —1 
on our parameter space. The results are shown in Figure |5] where 
the reconstructed EOS using the SNLS tracks the cosmological 
constant value at all redshifts, with 1— tail la errors of order 0.2 
at all z values. Because of the reduced freedom in w, the recon- 
struction collapses to the lower limit of the allowed w range, even 
after marginalization over the model values. In the case of the 
HST/GOODS data, a gentle rise of ■w{z) away from —1 is again 
observed. The larger error bars suggest that the entropy becomes 
important and that the value w ^ —0.8 to which the reconstruc- 
tion tends at redshifts z > 0.5 is mediated by the the mean value of 
the prior model m. As before, the best fit points remain very close 
to w — —1 at all redshifts. 

We can increase the amount of prior information by consid- 
ering the case where a constant prior model value m — —1 is 
used, see Figure |6] This is helpful in assessing whether the struc- 
ture observed in the HST/GOODS sample is strong enough to over- 
ride our entropic prior. The reconstruction from the SNLS data re- 
mains close to w = — 1 with errorbars of order 20% at all redshift. 
One has however to keep in mind that the tightness of the errors 
is partially helped by the supplementary information provided by 
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Figure 5. Reconstructed EOS (marginalized errorbars encompassing 68% 
of posterior probability) from SNLS data (top panel) and HST/GOODS data 
(bottom panel), marginalizing over a constant prior model but restricting 
the w range to — 1 ^ w ^ 0. The horizontal lines indicate the upper and 
lower bounds of the allowed w range (solid) in order to guide the eye. We 
find no significant deviation from ui = — 1 for the SNLS data set. For the 
HST/GOODS sample the reconstruction tends to drift to larger values at 
higher redshift. 



the entropic prior. This demonstrates how the use of cosmologi- 
cal constant as the model can be problematic as one can not say 
conclusively whether this indicates that the data is very strong or 
alternatively overridden by the entropy if no significant deviations 
from the model are observed. The result from the HST/GOODS 
data shows again a high value of w{z) being favoured in the third 
bin and subsequent collapse towards to the model m at higher red- 
shifts. The persistence of a deviation towards w > —1 at redshift 
z ~ 0.5 in the presence of the strong prior favouring the cosmo- 
logical constant suggests it is a real feature of the data. 

Finally, we also investigated the stability of the reconstruction 
against a change in the number of reconstructed components. Since 
the MaxEnt technique is designed to automatically deal with the 
structure in the data by adjusting the degree of smoothness of the 
reconstruction, we do not expect that a change in the number of 
bins would make a large difference in the reconstructed EOS. This 
is demonstrated in Figure|7] where the reconstruction analogous to 
the bottom panel of Figure [5] has been performed by halving the 
number of bins to A'^ = 5, without appreciable differences in the 
end result. 

Lastly, we investigate the dependency of our results on the 
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Figure 6. Reconstructed EOS from SNLS data (top panel) and 
HST/GOODS data (bottom panel), assuming an entropic prior model m = 
— 1 and a w range — 2 sC w ^ 0. The horizontal lines indicate the upper 
bound of the allowed w range (black) and the model m (red). The SNLS 
data are compatible with the model and show en'orbars of order 20% at all 
redshifts. The slight bump al z ~ 0.5 for the HST/GOODS data survives 
the entropic prior. 



0.5 



0.0 



-0.5 



-1.0 



-1.5 



-2.0 



1 1 ■ ■ 1 ■ ■ 1 ■ 1 1 ■ ■ 1 1 ■ 1 1 ■ 
HST/GOODS SNIa -I- other data 


■ A I 




; I 


I \ 


i 


T 




)K DIstrl 
ABest- 

1 


Jutlon Meon _ 
fit point 



0.0 



0.5 



1.0 
z 



1.5 



2.0 



Figure 7. Reconstructed EOS from the HST/GOODS data with 5 compo- 
nents in the reconstructed w. Compare with the bottom panel of Figure[3] 



model one chooses for the entropic prior. Up until this point we 
have assumed the prior model m to be constant with redshift such 
that it introduces a suitable degree of smoothing of time-dependent 
noisy features in the data. This choice reflects a specific belief in 
the true form of the EOS and given the large range of dark energy 
models on the market, it is important to assess the impact of a dif- 
ferent prior model m. Another popular class of models is given by 
an equation of state that is a smooth varying function of redshift, 
such as 



(1 



(34) 



Here the assumption is that the true EOS is a function of time and 
evolves sufficiently slowly such that it may be effectively charac- 
teri zed in a phenomenological way by the two parameters wq and 
u)i J Sahni & Starobinskrvl 2006h . This particular function is a good 
approximation to many dark energy models but clearl y it is limited 
to how well it can cop e with a rapidly evolving EOS jLiddle et aU 
I2OO6I) . Following e.g. Ilchikawa & Takahashil ( I 20061) , we impose 
the further constraint that the early Universe is matter-dominated, 
i.e. we impose the condition wq -\- w\ < on the prior model. 
In order to be as conservative as possible, we again marginalize 
over both prior model parameters, too and wi, as follows. The ex- 
pression for the entropy is again given by Eq. i22l . but the model 
coefficients are now given by 



a = (Co"' - 2) -I- (cr ~ 2) 



Zi 



+ 2. 



(35) 



Here (7™ and (7™ are the coefficient representing the prior 
model parameters. They are included in the hypothesis space, then 
marginalized over. 

We display our results using this slowly-evolving prior model 
in Figure [8] (with the model parameters marginalized over). This 
is to be contrasted with the analogous case of Figure [5] where 
the entropic model m is constant in redshift space. At low red- 
shifts the reconstructions for different prior models do not differ 
appreciably, with the most notable difference being a shift towards 
slightly higher values of w for the case of the SNLS data set. For 
the GOODS/HST programme data set, the peak in the third bin 
is recovered but is lower than in the previous constant m case. 
Again, the reconstructed EOS prefers higher values in the last few 
bins. At high redshifts, due to the accumulation of error, we ex- 
pect more entropic domination and hence a stronger weight of the 
prior model choice. However in this region dark energy becomes 
progressively less important and the ability of data to constrain the 
time-depedence of the EOS thus degrades considerably. It is en- 
couraging however that the low redshift behaviour of the recon- 
structed EOS agrees well for both marginalized prior models. This 
further strengthens our conclusions regarding the power and robust- 
ness of our technique. 



5 CONCLUSIONS 

Given that the dark energy models on the market are predominantly 
phenomenological, a reconstruction technique that does not require 
specifying a parametric form for w{z) would be an advantage. We 
presented such a method based on maximum entropy to reconstruct 
the equation of state of dark energy within a Bayesian framework. 
The principle of maximum entropy is invoked when assigning the 
Bayesian prior. This means that in the absence of genuine signal, 
the model that is most smooth, or that maximizes the infor- 
mation entropy is favoured, with the extent of this bias being de- 
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Figure 8. Reconstructed EOS (marginalized errorbars encompassing 
68% of posterior probability) from the SNLS data (top panel) and the 
HST/GOODS data (bottom panel) marginalizing over a prior model of the 
type {wq ,wi). The horizontal lines indicate the upper bound of the allowed 
w range (solid) and position of to = — 1 (dotted) in order to guide the eye. 
Comparison with Figure [5] where a constant prior model was used, indi- 
cates that the results from both data sets are robust against the choice of 
the prior model while exhibiting a mild dependence on this choice at higher 
redshift, where the constraining power of the data degrades and entropic 
domination sets in. 



termined by a regularizing constant that is automatically adjusted 
to the data. In our analysis we decompose into a sum of 

weighted orthogonal step-functions to facilitate the reconstruction 
of sharp features (provided binning is sufficient). Extensions of our 
analysis could easily include using alterative expansion functions to 
see whether other properties of the time-evolution of dark energy 
are detected or constrained. 

We find that the reconstruction of a dynamical w{z) using ar- 
tificial datasets of H{z) and Da{z) is very promising at low red- 
shifts but suffers from a bias towards the chosen default model at 
higher z. To combat this effect, the prior model was incorporated 
into our hypothesis space allowing the reconstruction of a model- 
independent distribution of w, with a manageable loss of accuracy. 
Once the technique was established and demonstrated, it was ap- 
plied to a combinations of the current cosmological datasets and 
two popular choices of prior models, namely a constant EOS and a 
mildly evolving w{z). 

Using a dataset including the current WMAP3 measurement 
of the CMB shift parameter, the baryon acoustic oscillation mea- 
surement and the HST Key project measurement of the Hubble pa- 



rameter in conjunction with the SN type la data from the SNLS 
project, we found that w — —1 in the redshift range ^ z ^ 1100, 
with errorbars depending on the prior and model assumptions. In 
the most optimistic case, where the data are supplemented by an 
entropic prior around w — —1, the error is of order 20% at all red- 
shifts. When the same dataset was instead supplemented by the SN 
sample from the HST/GOODS program, the results agree at low 
(z < 0.3) redshift. We found however that the reconstruction tends 
to prefer a value w > — 1 around z 0.5 with a significance be- 
tween 1 and 2(7, depending on assumptions. This shows the dangers 
of fitting a parametric form of 'w{z) to the data, in which case one 
is bound to miss possibly significant features in the measurements. 
The high-redshift behaviour of the EOS becomes increasingly dom- 
inated by the entropic prior and thus exhibits a mild dependence 
on the choice of prior model. We have investigated the correlation 
properties of our reconstruction, and identified a moderate anti- 
correlation among the first few redshift bins of our reconstructed 
points. 

The MaxEnt technique presented here improves on other 
methods designed to minimize noise artifacts in that the amount of 
information taken from the data is not determined by the analyst but 
rather dictated by the data themselves. The presence of real struc- 
ture rather than noise-induced complexity is indicated by the size 
of the error bars. The entropic prior adjusts the error bars when the 
information provided by the data is unreliable. Secondly, in the ab- 
sence of real information the reconstruction tends toward our most 
intuitive estimate of w{z) with suitably large variance. In conclu- 
sion the merits of this technique are that it is self-regulating in the 
sense that it allows the data to determine the amount of structure 
that is included. More importantly it does not require an inherent 
assumption of the functional form of the true equation of state. 

We hope that this technique will prove useful in deriving even 
stronger, model-independent constraints on the dark energy history 
from future, high-quality data. 
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